Local temperature in quantum thermal states 
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We consider blocks of quantum spins in a chain at thermal equilibrium, focusing on their properties 
from a thermodynamical perspective. In a classical system the temperature behaves as an intensive 
magnitude, above a certain block size, regardless the actual value of the temperature itself. However, 
a deviation from this behavior is expected in quantum systems. In particular, we see that under 
some conditions the description of the blocks as thermal states with the same global temperature 
as the whole chain fails. We analyze this issue by employing the quantum fidelity as a figure of 
merit, singling out in detail the departure from the classical behavior. As it may be expected, 
we see that quantum features are more prominent at low temperatures and are affected by the 
presence of zero-temperature quantum phase transitions. Interestingly, we show that the blocks can 
be considered indeed as thermal states with a high fidelity, provided an effective local temperature 
is properly identified. Such a result may originate from typical properties of reduced sub-systems 
of energy-constrained Hilbert spaces. Finally, the relation between local and global temperature is 
analyzed as a function of the size of the blocks and the system parameters. 



I. INTRODUCTION 

Since the early days of quantum theory, there has been 
a considerable effort to formulate the principles of ther- 
modynamics from a quantum perspective Recently, 
this approach has allowed to discover novel features of 
the nature of many-body systems, see for instance Refs. 
0, H, 0, [H, H|. In particular, careful analysis of the 
application of thermodynamics concepts to microscopic 
systems have shown the appearance of purely quantum 
features. In such systems, thermodynamic magnitudes 
may indeed loose their classical properties, giving rise to 
peculiar dependencies on some system parameters, such 
as size or total energy. A striking manifestation of this 
is given by the temperature, a magnitude considered to 
be intensive in classical thermodynamics. In a quantum 
scenario the temperature may not be well defined, re- 
sulting in the fact that subparts of thermal states may 
no longer be described as thermal states with the same 
global temperature as the whole system. More precisely, 
recent studies suggest that the intensive nature of the 
temperature may be lost not only in dependence of the 
size of the system subparts, as in a classical scenario, but 
also changing the temperature of the global system [j| H[ . 
It is then a relevant question to understand under which 
conditions the concept of temperature offers a correct de- 
scription of sub-parts of quantum thermal states. 

Although the main motivation of this work comes from 
a fundamental point of view, studying the limits of va- 
lidity of thermodynamics concepts may be relevant from 
a practical perspective. In fact, recent experimental 
progress in nano-sciences allows to access thermodynam- 
ical quantities, like temperature, at scales in which devi- 
ations from classical thermodynamics may become rele- 
vant 0]. For example, in Refs. [1,0 it has been pointed 
out that the breakdown of the concept of temperature 
might have consequences on thermometry and might be 
observed in experiments with spin chain compounds. 



In Refs. [1, 0] a set of conditions were established in 
order to assure that a large thermal system can be ap- 
proximated by a set of factorized blocks, each of them 
described in turn by a thermal state at the same tem- 
perature as the global system. Clearly, when such an 
approximation is valid, local measurements performed 
on a block provide results compatible with the global 
temperature. In other words, we do meet a situation in 
which temperature is intensive. However, local measure- 
ments on a single block do not provide any information 
about the correlations with the other blocks, which are 
simply disregarded (traced out) in the measuring pro- 
cess. Therefore, the approximation in Refs. 0,0] may be 
further relaxed still retrieving situations compatible with 
the concept of intensive temperature. As a matter of fact, 
one may only require that the block actually measured 
should be in a thermal state at the same temperature as 
the global one, not necessarily factorized with the rest of 
the system. This is the approach adopted in this work. 

The aim of the present work is to consider these is- 
sues exploiting ideas and tools recently emerged from 
quantum information science. This approach has the ad- 
vantage of giving a thorough description of the quantum 
systems under consideration, allowing in turn to charac- 
terize in detail the departure from the classical behav- 
ior. We consider a chain composed by n spins (where 
n can be taken in the macroscopic limit) at tempera- 
ture T and focus on the thermodynamical properties of 
a block composed by m spins (m < n) — i.e., the re- 
duced state obtained after tracing out n — m spins. In a 
standard thermodynamic setting m and n are taken large 
enough such that the interactions between the block and 
the rest of the system may be disregarded. As a con- 
sequence, the block can be well described by a thermal 
state at the same temperature as the global one. How- 
ever, such a picture may brake down for blocks of small 
size and strong interactions. As said, whereas for clas- 
sical systems this breakdown has no dependence on the 
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temperaturejin quantum systems a temperature depen- 
dence arises Q . This is the scenario that we will consider 
here. To face this problem we use the quantum fidelity 
as a figure of merit. The latter quantifies the amount of 
statistical distinguishability between two quantum states 
and its properties have been widely studied in quantum 
information science 10]. Recently, the fidelity has been 
used both as an indicator of phase transitions in spins 
systems [ll], [IH , and to study the emergence of thermal 
states in contiguous blocks of some spin-system ground 
states 13]. As shown in these works, the fidelity has 
revealed to be a particularly sensitive figure of merit in 
this framework. Here, we use it to check whether the re- 
duced states of the spin chain under consideration can be 
well approximated by an m-particlc thermal state. The 
sensitivity of this approach will be assessed by investigat- 
ing regions in the phase space near the zero-temperature 
critical points. 

We also relate our investigations to recent findings in 
the foundation of statistical mechanics - again inspired 
by quantum information concepts. In Refs. 0, Q the 
authors address typical properties of the reduced states 
obtained by tracing out a huge amount of degrees of free- 
dom from a constrained pure system. Building on previ- 
ous results on properties of quantum states in large di- 
mensional systems [l4| , it was shown in Ref . @ that the 
reduced state of a big system (including environment) 
satisfying an operator constraint is basically the same 
for almost any pure state of the system. In the particu- 
lar case in which (i) the operator constraint is related to 
the energy of the whole system and (ii) the interaction 
between system and environment is small, this typical 
state can be shown to correspond to a canonical thermal 
state @, Q , as already pointed out in the early days of 
quantum mechanics [151 ]. 

Motivated by these results we study the reduced states 
of the chain from a thermodynamical perspective, de- 
scribing it with only a few physical magnitudes, an effec- 
tive temperature in our case. In other words, we check 
whether the reduced states of thermal systems maintain 
some sort of canonical typicality. Recall that the results 
in Refs. [5j, [6| show that a canonical thermal state for 
the reduced system is obtained when the interaction en- 
ergy between the considered parts is negligible. As said, 
here we will consider settings in which such condition 
is not fulfilled, hence the validity of these typicality re- 
sults is by no means guaranteed. Specifically, we pro- 
ceed by identifying the thermal state of m particles that, 
subject to the local interaction inherited by the whole 
Hamiltonian, is closer to the actual reduced state of the 
whole chain. We thus define an effective local temper- 
ature (the only free parameter to adjust) that can be 
compared with the global temperature of the whole sys- 
tem. Depending on the parameters of the Hamiltonian 
and the subsystem size we find situations where the local 
temperature is no longer equal to that of the global sys- 
tem, in accordance with the above mentioned analysis. 
However, the description of the reduced states as ther- 




FIG. 1: (Upper panel) Reduced m-particle state p m (/3) of a 
chain composed by n particles in the thermal state Q„(f3) at 
temperature T — 1//3. The state Sl n (/3) is built from an n- 
particle Hamiltonian H n . (Lower panel) Thermal state Q m (/3) 
of m particles for the same local interaction H m but referred 
to m particles. 

mal states is an extremely good approximation, as shown 
by the corresponding fidelity. This is remarkable since, 
as said, we are not in the standard conditions considered 
in Refs. [E,@]- Our results then suggest that some form of 
canonical typicality may still be present, even if the inter- 
action between system and environment is not negligible. 
In particular, we find that below some threshold temper- 
ature the local temperature may become higher than the 
global one, and that the reduced states can have some 
finite temperature even when the global system is in the 
ground state. Furthermore, as one may expect, the local 
and global temperatures tend to coincide increasing the 
size of the subsystem. 

The structure of the paper is the following. After in- 
troducing the systems under consideration, we calculate 
in Section [TT] the fidelity between two-particle blocks and 
the corresponding thermal state. This allows to study 
the conditions under which the temperature ceases to be 
an intensive magnitude. Then, in Section IIII1 we derive 
the local effective temperature of the two-particle blocks, 
as well as the fidelity of this description. The dependence 
of the local temperature on the size of the blocks is stud- 
ied in Section ITVl where blocks consisting of more than 
two spins are studied with the help of the Matrix Prod- 
uct States formalism. Finally, we summarize our results 
in Section fVl and point out some conclusive remarks. 



II. INTENSIVE TEMPERATURE ANALYSIS 

We start by assessing the applicability of the concept 
of intensive temperature. As said, we consider a large 
spin system at temperature T — 1//3 and focus on a sub- 
block of spins. We then check whether the latter can 
be well approximated by the thermal state at tempera- 
ture T given by the interaction inherited from the whole 
system Hamiltonian. By analyzing the relation between 
the reduced states of thermal systems and thermal states 
themselves, we assess the validity of such estimation. 

Consider the setting depicted in Fig. [T] On the one 
hand, we construct the canonical state of n particles 
f2 n (/3) = exp (— PH n )/Z n , where H n is a local interact- 
ing Hamiltonian and Z n is the partition function. In 
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order to check the thermal properties of the subparts of 
this state we trace out a part of it, obtaining the state 
of m particles, p m {P) = Tr„_ m [fi„(/3)] (throughout this 
work we will focus on systems composed by a number 
of particles n much larger than m). On the other hand, 
we directly construct the canonical state of m particles 
^m(/3) = exp(—f3H m )/Z m . We then compare these two 
density matrices using the fidelity measure [l6| 

1/2 



F[n m (0),p m (/3)] = Tr ( Vn m {j3)p m {j3)Vn m (P) 

(1) 

Throughout this Section both fl n (/3) and O m (/3) are set 
to the same temperature, since we want to identify when 
the temperature is intensive. 

We apply the above considerations to a 1-D spin chain 
characterized by the Ising Hamiltonian in a transverse 
field with open boundary conditions 
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where <Ji are the Pauli matrices, and h gives the strength 
of an external magnetic field. The system experiences a 
quantum - i.e^, zero temperature - phase transition when 
h = h c = 1 1171 . The two-spin correlation functions are 
given by 18l.ll9|| 
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The parameter r sets the distance between the particles, 
e.g. r = 1 means two neighbor particles. These correla- 
tors are calculated for chains in the thermodynamic limit 
(i.e., n — * 00). Using these formulae, one can compute 
the reduced density matrix for the 2-spin system f>2{/3), 



(8) 



without the explicit construction of the global thermal 
state of n particles . 



The minimal size for which both the terms in the 
Hamiltonian @ contribute to the construction of VL m (/3) 
is for m = 2. We then devote much attention to this 
first non-trivial case of two-spin blocks. Furthermore, it 
is reasonable to expect that, if for such small blocks the 
temperature is intensive, it will be intensive for larger 
blocks as well. This intuition will be confirmed in Sec- 
tion [iVl where the size dependence of our considerations 
will be analyzed. 

Before proceeding with the results, let us comment 
about our choice of the local Hamiltonian H2, from which 
the reference state f22(/3) is derived. With this choice we 
are at first sight disregarding the interaction between the 
2-spin block and the rest of the system. Whereas this 
may be easily justified in the case of large m and n, it 
deserves to be clarified in our setting. In fact, we consid- 
ered various strategies in order to take into account the 
interactions at the border of the block, in particular a 
mean-field approach, and the limit of high temperatures. 

First, one may add to H2 a correction term which takes 
into account the surrounding spins of the block following 
a mean-field approach. The idea is to replace the opera- 
tors at the boundaries by their mean value (<t x ), so the 
corresponding effective two-particle Hamiltonian reads 



H'2 — 7: a x 



u 2 



(9) 



The boundary correction term however turns out be zero 
for any finite temperature, as (a x ) = 0. This can be seen 
by considering that the canonical distribution fi n (/3) in- 
herits the symmetries of the Ising Hamiltonian, in par- 
ticular the global spin-flip symmetry U = As a 
consequence, one has that [U, f2 n (/3)] = 0, which implies 
that 



„PiG8)] = 



(10) 



for any finite (3. Considering now that pi(f3) can be ex- 
panded in terms of the Pauli matrices, Eq. (|10l) imposes 
that (<j x ) = for any T > (see also Ref. [19| ). How- 
ever, recall that a symmetry breaking can occur at T = 0, 
resulting in (a x ) 7^ 0. 

Second, one could consider a correction valid for high 
temperatures. A first order expansion for — * imme- 
diately reveals that in this limit P2(P) is given by 



Tr„_ 2 ( — - pH n 



-^2 0H% 



(11) 



since the Pauli matrices appearing in the interacting 
terms of H n are traceless. The expression above for foiP) 
coincides with the high-T expansion of SI2 (/?) . In other 
words, for high temperatures the standard situation is re- 
trieved, and no correction has to be taken into account. 

In order to further clarify our choice of the local Hamil- 
tonian, let us consider the classical antiferromagnetic 
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one-dimensional Ising model, given by: 

n-i 

/#=5>V+ 1 , (12) 

»=i 

where s l — ±1. Notice that the model above gives the 
classical limit of H n for the case h = [13]. It is worth- 
while to briefly consider this example since it gives a rele- 
vant classical model whose local Hamiltonian contains no 
correction due to the boundary terms. To show this, let 
us consider the thermal state associated to Hamiltonian 
([12"]). described by a thermal probability distribution: 

P( S 1 ,.., S ™) = ^exp[-/3^ 1 ] , (13) 

where Z n is the partition function. A straightforward 
calculation shows that, by summing over all the possible 
configuration of the n — 2 spins surrounding an arbitrary 
2 spin block, one obtains the following distribution func- 
tion: 

^2 

where k is arbitrary (1 < k < n), — s k s k+1 is the lo- 
cal Hamiltonian and Z2 its respective partition function. 
As said, we see that for this kind of classical system no 
correction due to the boundary terms is expected. Fur- 
thermore, notice that Eq. ([14")) implies that the tempera- 
ture is intensive. However we should stress here that this 
result is not true for a generic classical system. In conclu- 
sion, given the considerations above, it seems reasonable 
to consider the Hamiltonian J7 2 ~~ without any additional 
correction - to build the reference state f^C^)- We will 
adopt this choice throughout the paper. Finally, let us 
notice that the same considerations apply for any size m 
of the block (see Sec. HVJ). 

The results for the fidelity F[n 2 (/3), h{0)\ are plotted 
in Fig. [21 We can see that it turns to be a non-trivial 
function of the external magnetic field h and the system 
global temperature T. In particular we observe a high 
fidelity above a certain temperature, for any value of h. 
This confirms the intuition that for high temperatures 
the classical behavior should be recovered. Namely, the 
reduced states are well approximated by thermal states at 
temperature T . On the other hand, as we lower the tem- 
perature, the fidelity can drop to values sensibly lower 
than one, thus indicating that the standard thermody- 
namic description of the reduced state is no longer ac- 
curate. We then recover the main feature of the results 
given in Q: the validity of the concept of intensive tem- 
perature depends on the temperature itself, a behavior 
with no classical analogue. Thanks to the sensitivity of 
the fidelity measure, we can also analyze in details such 
a behavior, as can be seen for low temperatures. In this 
case, the fidelity turns to be equal to one for h » 1 and 
h <C 1, as can be expected recalling that in both cases 
the ground state is factorized (all the spins are aligned 




FIG. 2: Fidelity F[tt 2 (f3), p 2 (/3)} between the states fi 2 (/3) 
(the thermal state of two spins at temperature T = 1/(3) and 
p2(/3) (the two-spin reduced state of a global thermal system 
at temperature T) as a function of (3 and the local magnetic 
field h of the Ising Hamiltonian ([2]) . The temperature ceases 
to be an intensive magnitude when the fidelity is significantly 
different from one. 



along the same direction) [17j. On the other hand, for 
intermediate magnetic fields, the fidelity drops to values 
sensibly lower than one, showing a minimum which de- 
pends on the temperature. We recall that the model in 
Eq. ||3J| has a critical point when h c = 1. 

Before proceeding further, we now explicitly assess the 
sensitivity of the fidelity in this scenario. For this pur- 
pose we numerically evaluated the first derivative of the 
fidelity with respect to h at fixed (3. We see in Fig. [3] 
that, as (3 increases, the maximum of the derivative gets 
more pronounced. This behavior reflects the fact that the 
changes in the ground state become sharper and sharper 
as the system approaches criticality. As a consequence, 
subparts of the system change sharply as well and the 
approximation of them with a thermal state becomes 
more sensitive to small changes of h. In turn, the ac- 
tual value of the maximum in Fig. [3] also increases as 
the temperature decreases. This behavior proceeds un- 
til we reach some /3, which corresponds to some effective 
ground state, and below which additional changes are 
hardly observed. This can be understood by looking at 
the temperature dependence in the correlators in Eq. ^ . 
For large /3, one can well approximate tanh(i/3o;0) » 1, 
so the reduced states become almost independent of the 
temperature. The derived magnitudes computed from 
the reduced states exhibit for this reason minor changes 
below some temperature. 

The above considerations can be extended as well to 
two spins separated by r particles in the chain. In partic- 
ular, we considered a non-contiguous block of two distant 
spins which we denoted by f>2.r{P)- We reconstruct the 
density matrix of such a block using again Eq. (jHJl and 
considering the explicit dependence on r of the correla- 
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FIG. 3: First derivative of the fidelity F[fi 2 (/3), fatf)] with 
respect to h and for fixed values of f3 (/3 — 50, 100, 500). The 
behavior turns singular at the critical point h c for values of /3 
closer to the ground state. 



tors [see Eq. ([BJ]. Clearly, it is no longer interesting now 
to compare p2,r with U.2 (/3) . A much reasonable strat- 
egy is instead to construct a thermal state f2 r +i(/3) com- 
posed by r + 1 spins and trace out all the particles but 
the two extremal ones. We denoted the two-spin state 
so obtained as fl2,r(P)- Then, we can compare fi2,r(/3) 
and f>2,r(l3) by calculating the fidelity F[0 2 ,r(/3), f>2,r(P)]- 
The results for different values of r are shown in Fig. [4] 
Clearly the fidelity is higher with respect to r = 1. How- 
ever, some qualitative features already observed in the 
case r = 1 are still present for larger r. In particular, the 
fidelity turns equal to one for /i> 1 and h <C 1, whereas 
for intermediate magnetic fields it shows a minimum. 

It is worth summarizing at this point the two main 
features disclosed by the above analysis. Namely, i) the 
intensive nature of the temperature depends on the global 
temperature of the system itself, and ii) the temperature 
ceases to be intensive in a limited region of the phase 
space, namely for intcrmadiate magnetic fields around 
the zero-T critical point. As we will mention later on, 
we obtained similar results also in systems different from 
the one given by Eq. @. 



III. EFFECTIVE LOCAL TEMPERATURE 

We have seen in the preceding section that, for some 
values of the Hamiltonian parameters and the tempera- 
ture, a two-particle block of a thermal state may be dif- 
ferent from a two-particle thermal state under the same 
Hamiltonian at the same temperature. Now, we focus 
on the reduced state p m and look for a valid descrip- 
tion of it in terms of only a few thermodynamic magni- 
tudes. A key point in the attempt to describe the reduced 
states is the fact that, under some natural circumstances, 
they will become with high probability close to a ther- 
mal state [5, 6] . This follows from the general structure 



FIG. 4: For fixed (3 = 10 we plot the fidelity 
F]p,2, r {0), p2,r{P)] between the states Q2,r(/3) [the reduced 
state of two particles distant r obtained from a thermal state 
Q, r +i(f3)] and p2, r (J3) [the reduced state reduced state of two 
particles distant r obtained from a thermal state Q n (f3) with 
n — > 00] for r = 3, 5, 7 (bottom to top). Recall that r denotes 
the distance between two spins. 



of the Hilbert space of the global system, and the re- 
striction over the total energy 2l| . In this direction, the 
results in Refs. 0, @ are quite general but their applica- 
tion to a specific system has to be taken cautiously. In 
particular, a thermal state is recovered when the interac- 
tion between the parts of the system under consideration 
is negligible with respect to the interactions inside the 
parts themselves. This condition, however, may not al- 
ways be satisfied. When the parts taken into account are 
composed of two spins only, as analyzed in the previous 
section, it is then non trivial that the reduced states are 
actually in a thermal state of their respective Hamilto- 
nian. Nonetheless, we will see that this is indeed the case 
in the majority of the circumstances. 

Following these considerations we check whether the 
reduced state is a thermal state as well, but at an ef- 
fective local temperature different from the one of the 
global system. That is, even if the temperature may not 
be intensive, the reduced states may still have a sim- 
ple thermodynamical description. To this end, we con- 
sider a generic m-particle thermal state fl m (f3') and look 
for the effective optimal temperature (3 for which such a 
state better describes the actual m-particle reduced state 
p m {[3). As outlined in the previous Section, we assume 
that the local interaction H m is the same as in the global 
system H n (actually H m is obtained from H n by tracing 
out the n — m disregarded spins), we only have to adjust 
the temperature. In order to check the quality of this 
description, we compute the fidelity between the reduced 
state and the reference thermal state, and then optimize 
the temperature of the latter. That is, once we compute 
the traced state p m ((3) for any given (3, we optimize over 
the parameter (3' the fidelity function F[f2 m (/3'), p m {(3)}. 
In this way, as said, we identify an effective local temper- 



FIG. 5: The optimized inverse temperature (local temper- 
ature) of a reduced subsystem of two spins as a function of /3 
(/3 = 1/T, where T is the global temperature of the system), 
for different values of h. 

ature (3' = (3 for the reduced state of the system. 

Operating similarly as in the preceding section, we 
study for m = 2 the relation between f3 and (3. As 
shown in Fig. [5j we have that j3 w (3 for low values of 
/3 (i.e., the local temperature is the same as the global 
one). Thus, this range of temperatures can be identified 
as a classical regime, where both the local and the global 
temperature coincide, and the temperature behaves as 
an intensive magnitude. By lowering the temperature 
this equality stops to hold, and the value of (3 saturates. 
For even lower temperatures of the global system, the 
reduced state keeps the same effective temperature, even 
in conditions of zero-T where the whole chain is in its 
ground state. As above, this is due to the fact that 
for sufficiently small T, the reduced states hardly change 
with temperature. Notice also that the relation /? j3 
holds even for small temperature when the magnetic field 
vanishes {i.e., h <C 1). This is true also for h ^> 1 (not 
shown in the picture) and can be understood in view of 
the considerations made in the previous Section. 

Performing the optimization over the local tempera- 
ture we improved dramatically the previous values of the 
fidelity, as can be seen comparing Figs. [21 and \§\ [221 ] . 
As a matter of fact, after the local optimization the fi- 
delity is everywhere almost one, attaining the minimum 
of F[fl20),p 2 (/3)} « 0.975 for h w 0.6. Such result im- 
plies that a local temperature (3 may be defined for almost 
every (3 and h, even if such local temperature is no longer 
intensive. In other words, the findings in Refs. [H, [|| can 
be applied, at least approximately, also in cases like the 
one studied here. This is in a sense surprising, since here 
the assumptions used to derived the results of Refs. 0] 
are no longer satisfied. In particular, the block we are 
considering is composed by only two particles, implying 
that the interaction between it and the rest of the system 
cannot be disregarded a priori. Let us notice moreover 
that the definition of a local temperature is not com- 



FIG. 6: Optimized fidelity F opt = F[Q 2 (/3), p 2 (/3)] between 
the thermal state Q2 ($) at an optimized inverse temperature 
(3 (local temperature) and the reduced state f>2{(3) of two spins 
at inverse temperature fj (global temperature). Compared 
with Fig. [2] the values of the fidelity are sensibly higher. 
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FIG. 7: Value of a ^ as a function of the local field h (/3 = 
100,250, 1000). 



pletely satisfactory for intermediate magnetic field h at 
low temperature. This resembles the results obtained in 
the previous Section even if, as said, the fidelity here ob- 
tained is much higher. We will commment on this in 
Sec. El 

In order to check the sensitivity of the local tempera- 
ture (3 to small changes of the parameters we studied the 

derivative as a function of the local field h which 

drives the quantum phase transition. We plot the re- 
sults in Fig. [71 where a singular value around the criti- 
cal point appears. This shows that minor changes in the 
global temperature will strongly change the local thermal 
properties of the reduced subsystem. The comparison of 
Fig. [3] and [7] shows that the influence at finite T of the 
quantum phase transition is more pronounced for local 
effective temperatures. 
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FIG. 8: Fidelity F[p 2 (f3), p 2 (/3+A/3)] between the two reduced 
states of a system at inverse temperature f3 and /3+ A/3 respec- 
tively (A/3 = 0.3, 0.5, 0.7 from top to bottom). The magnetic 
field is set to h — 0.1. 

We complete the study of two-particle reduced states 
comparing two of them at slightly different temperatures. 
In particular, we compute the fidelity F[p 2 (P), p 2 {P + 
A/?)] between two reduced states at temperatures /3 and 
f3 + A/3, as a function of the global system temperature 
(3. Here p 2 {(3) is again a two-particle reduced state of a 
chain at global temperature (3. In Fig. [8] we show the 
fidelity as a function of (3 for different values of A/3. The 
value of (3 for which F[p 2 ((3), p 2 (f3 + A/3)] = 1, combined 
with the preceding results, suggest that above a given (3 
the states are equal and independent of the temperature. 
This is a direct indication of the asymptotic character 
of the reduced states p 2 (f3) for low temperatures, as al- 
ready expected from Eq. ([6]). In addition, the values of 
the fidelity for different A/3 appear converging to one at 
a value (3 consistent with the saturation of the local tem- 
perature found in this Section. 



IV. ANALYSIS FOR INCREASING BLOCK SIZE 

In the previous Sections we mainly discussed reduced 
blocks consisting of two spins, not necessarily contiguous. 
Now, we extend our considerations to larger blocks. Re- 
call that in the case of two-particle blocks we considered 
the whole system in the thermodynamic limit (n — ► oo) 
and were able to construct the density matrix of the block 
via the two-body correlators in Eqs. ©-([5]). Even if such 
a procedure can be in principle extended to higher or- 
der correlators (for the 3-spin correlators see for instance 
[231), it is more convenient to use a different strategy 
when blocks of arbitrary length are considered. Specif- 
ically, we use the formalism of Matrix Product States 
(MPS), a numerical tool which has been shown to de- 
scribe with high precision ground and thermal states of 
ID local Hamiltonians. Using the MPS formalism for 
mixed states [2J, we construct thermal states and 



compute correlation functions in an efficient way. Within 
the adopted numerical approach, we cannot work directly 
in the thermodynamical limit. However, the MPS formal- 
ism allows to simulate systems large enough to accurately 
reproduce the thermodynamical limit. We checked that 
systems composed by around 50 spins already suffice to 
obtain results indistinguishable to the corresponding ones 
in the thermodynamical limit. We will explicitly show in 
the following such a comparison. 

The interest in considering blocks of increasing size 
is motivated by the following observations. On the one 
hand, when the block increases in size the interaction 
between the latter and the rest of the system becomes 
less relevant when compared with the interactions in- 
side the block itself. This may intuitively lead to guess 
that the rest of the system plays a negligible role in the 
thermalization process. In other words, each block in 
which the global system may be divided could be consid- 
ered as a system on its own, thus thermalizing with the 
global environment independently from the rest. In such 
a scenario the temperature would clearly be intensive, for 
large enough block size. On the other hand, there are sit- 
uations in which it is not possible to disregard any part 
of the whole system, in particular when the latter is in 
a highly correlated state (e.g., near a phase transition). 
Thus, an analysis of the intensive nature of temperature 
as a function of the system size can identify which of 
these two tendencies prevails and in which setting. 

Before proceeding, let us briefly recall the MPS for- 
malism (see [24| for more details). This representation is 
based on a set of matrices A s k of size D x D used to write 
a quantum state as 

d 

\*mps)= £ Tr(A?...A£)\s 1 ...s n ). (15) 

Sl ...s n = l 

We can write mixed states in the MPS formalism intro- 
ducing at each position an ancillary system of dimen- 
sion d. Thus, the thermal state p is written as a pure 
state in a Hilbert space of larger dimension as 

d d n 

\P)= E E Tr(Y[A s k ^)\ Sl ... Sn )\ ai ...a n ). 

si . ..s n a\...a n k—1 

Using this purification one can recover the thermal state 
tracing out the ancillary systems a^, p = Tr a (|p)(p|). 

With this representation we can compute expectation 
values efficiently using the relation 

(Or . . . O n ) p = Tr(E h0l ■ ■ ■ £„,oJ , (17) 

where Ek,o = J2s s> ( s 'l<^l s )-^fc' S > an d the set of matrices 
Mfc are the result of tracing the ancilla states 

d 
a=l 
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FIG. 9: The local inverse temperature (3 of a reduced subsys- 
tem of m spins as a function of the global inverse temperature 
of the system (the magnetic field is set to h = 0.8; the ma- 
trices size for the MPS algorithm is set to D = 15). From 
bottom to top we set m = 2,3, 4, 5, 6. The dashed-dotted line 
shows = for reference. We see that for larger blocks the 
local temperature gets closer to the global one. The actual 
region in which the temperature is intensive does not seem to 
be affected by the size of the block (m > 2). 

With an initial MPS representation of a state at (3 = 0, 
which corresponds to the completely mixed state 1, we 
obtain the state p cx exp(—0H) using 

& - m = (e -A t ff ) M 1(e -A t ff ) M j (1Q) 

where At = 0/2M. In this way the thermal state at 
temperature is the result of the evolution in imaginary 
time of the completely mixed state, an evolution that can 
also be performed in an efficient way by means of a Trot- 
ter decomposition of the time evolution operator. Using 
this numerical technique we can extend the calculation 
of the correlation functions up to m = 6, the computa- 
tional limit being the construction of the state p from the 
corresponding correlators as in Eq. (jHJ). 

We plot in Fig. |S] the effective local temperature, 0, as 
a function of global one, (3, for reduced systems of size 
to = 2 — 6. The first feature that one may note is that 
the local optimized temperature tends to get closer to the 
global one as the block size increases. This is in agree- 
ment with the trivial fact that for large block size the 
local and global temperature coincide, as the block itself 
coincides with the whole system. Furthermore, as men- 
tioned above, as the block size increases the interaction 
between the latter and the rest of the system becomes 
less relevant. Thus, the block will be in a thermal state 
with high probability, as a consequence of the results ex- 
posed in Refs. [!, [||. For each curve depicted in Fig. 
we checked explicitly the fidelity between a thermal state 
at temperature (3 and the actual reduced state. We ob- 
tained that, all along these curves, the fidelity is very 
close to one, confirming that the results in Refs. [H, Q 
can be applied to the finite systems considered here. In 



FIG. 10: We plot the local inverse temperature (3 reached 
by the reduced subsystem of m spins as a function of the 
magnetic field h (the global inverse temperature is set to (3 = 
15). From bottom to top the size of the block increases (m = 
2,3,4,5,6). 

turn, this points to the fact that the systems are already 
near to the thermodynamic limit, at least for what con- 
cerns the properties of their reduced states. However, the 
fact that the blocks are thermal does not mean that the 
temperature is intensive. As said, Fig. [9] shows that the 
local temperature get closer to the global one for larger 
blocks, nevertheless the actual region in which the tem- 
perature is intensive does not seem to be much affected 
by the size of the blocks. This suggests that the quan- 
tum fidelity analysis reveals fine features of the systems 
that cannot be understood by considerations regarding 
only the energy balance between the various parts of the 
system. In particular, as mentioned above, correlations 
should be considered in order to clarify such a behavior. 

We also notice that for each value of to a different sat- 
uration value of is obtained. The latter corresponds to 
the effective temperature of the reduced states when the 
whole system is at zero T. Thus for the ground states of 
the considered systems, the temperature is clearly not in- 
tensive, as expected. On the other hand, for high enough 
temperature, we can identify again a classical regime — 
i.e., a regime in which the temperature is intensive. In 
general, outside this classical region we observe that the 
local temperature is sensibly higher than the global one. 

Before concluding this Section let us consider the influ- 
ence of the zero-T properties of the system at small tem- 
peratures. As seen in the previous Section (see Fig. [7]), 
we expect that at small temperatures the features of the 
system ground state can be revealed by an analysis of 
the local temperature as a function of the local field h. 
The results for = 15 are plotted in Fig. [TUJ In accor- 
dance with the results reported in Fig. [51 as the block size 
increases the local temperature gets closer to the global 
ones for any value of h. Furthermore, notice that, as the 
block size increases, the minimum of (that identifies 
when the temperature deviates mostly from an intensive 
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behavior) get closer to h c . This behavior around the criti- 
cal value of the field h = 1 suggests that zero-T properties 
of the system ground state have an important influence 
on the thermal state of the block. 



V. CONCLUSIONS 

With the study of many-body systems using thermo- 
dynamic quantities, one may obtain global properties of 
a system by measurements performed only on a local (re- 
duced) part of it. However, depending on the conditions 
of the physical system under consideration, some of the 
fundamental assumptions of statistical equilibrium may 
fail and a valid thermodynamical description cannot be 
possible. We have identified this situation in spin chains 
at low temperatures, where the reduced states are not 
thermal states at the same temperature as the global 
system. In particular, we have seen that the tempera- 
ture ceases to be an intensive magnitude in dependence of 
the global temperature itself, a feature without classical 
analogue. This result is in accordance with the findings 
reported in Ref. [jj. Notice, however, that here we made 
no assumption on the correlations between the blocks in 
which the global system may be decomposed. In fact, 
we focused only on the reduced local state of the sys- 
tem. Our approach is motivated by an operational view- 
point, since an actual measure performed on a local part 
of a system gives no information about correlations with 
other parts of it. 

Remarkably, we have seen that an effective thermo- 
dynamical description of the reduced states is often still 
possible, despite the temperature may not be an intensive 
magnitude. In particular, the resulting reduced states 
can be described with high accuracy as thermal states at 
an optimized temperature. Interestingly, this is valid i) 
without any additional assumption on the interactions, 
and ii) even if we considered blocks of particularly small 
size (i.e., composed of 2 to 6 spins). This local temper- 
ature becomes of course equal to the global one in the 
regime of high temperatures, where we recover the clas- 
sical behavior. However, by lowering the temperature 
of the global system, the local temperature saturates at 
some point. In other words, depending on the system 
parameters, the local temperature may be different from 
zero even when the global system is in the ground state. 
The dependence of the local temperature on the size of 
the blocks has been studied too. As one may expect, the 
local temperature get closer to the global one for larger 
blocks. 

We pointed out that the departure from the classical 
behavior is more pronounced at low temperatures and for 



intermediate values of the local magnetic field, around 
the zero-T critical point (this is true concerning both the 
intensive behavior of temperature and the definition of 
a local temperature, as can be seen in Figs. [5J 2J and 
[6J. Notice that this is the region of parameters where 
quantum correlations are expected to be stronger. For 
example, at zero temperature the pair- wise entanglement 
between two spins shows its higher values in this region 
[l9l |2H | , as well as the entanglement between a block of 
spins and the rest [27| . Though the existence of a fi- 
nite block entropy may suggest that entanglement plays 
a crucial role in all these results, preliminary studies indi- 
cate that the relation between quantum correlations and 
the definition of a local temperature is nontrivial. Thus 
it will be a subject of further studies the way in which 
quantum and classical correlations lead to the bahaviours 
identified here. 

Let us mention here that we have analyzed also other 
Hamiltonian systems. In this paper we have shown a 
study of the transverse Ising model, but the extension 
to a generic XY interaction leads to similar results and 
conclusions. Applying the MPS techniques we can ex- 
tend the analysis to spin models with arbitrary inter- 
action, such as spin chains with Heisenberg interaction. 
Furthermore, we analyzed chains consisting of harmonic 
oscillators with quadratic interaction (harmonic chains). 
In all these cases, the obtained results were qualitatively 
very similar to those shown here. 

In order to assess the canonical character of the re- 
duced states we employed the quantum fidelity, a quan- 
tity extensively used in quantum information science. We 
have seen that the fidelity is particularly suitable when 
studying the limits of applicability of the concept of in- 
tensive temperature. In particular, harnessing the sensi- 
tivity of quantum fidelity, we have observed how the low 
temperature features of the systems influence our results. 
This suggests that the quantum fidelity analysis reveals 
fine features of the systems that cannot be understood 
by considerations regarding only the energy balance be- 
tween the various parts of the system. In particular, as 
said, future investigations should consider the role of cor- 
relations in order to clarify such a behavior. 
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